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Abstract. We consider a possibility to unify the methods of regularization, such as the 
^s ' renormalization group method, stochastic quantization etc., by the extension of the standard 

t^|- ■ field theory of the square-integrable functions 4>{b) € L 2 (M' i ) to the theory of functions that 

CN \ depend on coordinate b and resolution a. In the simplest case such field theory turns out 

to be a theory of fields 4> a {b, ■) defined on the affine group G : x' — ax + 6, a > 0, x, b £ 

R d , which consists of dilations and translation of Euclidean space. The fields </> a (6, •) are 
(^-») . constructed using the continuous wavelet transform. The parameters of the theory can 

[ — ' explicitly depend on the resolution a. The proper choice of the scale dependence g = g(a) 

makes such theory free of divergences by construction. 
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23 ■ 1 Introduction 

In many problems of field theoretic description of infinite-dimensional systems the continuous 
y\ ' description of fields and propagators faces the problem of infinities in loop integrals. The 

underlying physics is that field theoretic description includes integration over all microscopic 
scales smaller than the size of the system. In the simplest case small-scale fluctuations can be 
just averaged to get the classical equations for macroscopic fields. The best known example 
is the laminar flow hydrodynamics. In more complicated cases interaction strength (and even 
type) may depend on the scale. The examples are: turbulent fluid flow, critical phenomena, 
elementary particles interaction, etc. 

For this reason we need a measure of integration that works better than Euclidean integration 
does. For some cases such measures are observed experimentally: these are fractal measure of 
energy dissipation in turbulent field P3EJ; fractal mass distribution of percolation clusters, etc. 
For high energies, when the fractal distribution of the field or order parameter can not be 
measured experimentally due to the uncertainty principle, the search for a better measure of 
integration remains a mathematical problem. Fortunately, all these phenomena share a common 
important feature - the self-similarity, or scaling. The percolation processes in nanoelectronics, 
the turbulent velocity field, and the nucleon scattering - all display the same picture if being 
observed at different scales. The self-similarity has initiated at least two important regularization 
methods for field theory models: 

• The renormalization group (RG) method that makes use of substitution of initial fields 
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4>(x) G L 2 (M ) by the scale-truncated fields 

<%)(*) = / e"**)^, (1) 

1 A > J\k\<\ (2vr) fl 

makes the coupling constants dependent on the cut-off momentum A, and requires that 
the final physical results should be independent of the introduced scale 

Kd\ (Physical quantities) = 0. 

• Use of random processes that are self-similar by construction for the regularization of 
quantum field theory (QFT) models. This is known as stochastic quantization. 

In present paper we summarize the both methods by introducing the fields <p a {b) that expli- 
citly depend on both position b and the scale (resolution) a. If 4>{x) is considered as the wave 
function of a quantum particle, then the normalization f \(p a (b)\ 2 dfi(a, b) = 1, where d/j,(a,b) 
is the appropriate Haar measure on affine group, expresses the fact that the probability of 
locating the particle anywhere in space — oo < b < oo changing the resolution < a < oo is 
exactly one. Thus, the RG symmetry related to the change of the scale o extends the symmetry 
of the theory by allowing the change of parameters with scale. Technically, construction of the 
scale-dependent fields (f> a (b) is performed by using continuous wavelet transform (CWT). 

It should be noted that since field theoretic calculations are usually performed with the help 
of the Fourier transform the cut-off momentum A corresponds to the minimal coordinate scale 
L = =£■. However, the Fourier transform is not localized in coordinate space. For this reason we 
denote the fields obtained from cj)(x) by truncation of the Fourier harmonics with \k\ > ^f- Q 
as 4>il)( x )'i i n contrast to it the localized view of a function <f> at the position x and the scale 
a is denoted as (p a (x). The mathematical meaning of the latter will be explained hereafter in 
terms of wavelet transform. The same thing concerns the parameters of field theory models - 
masses, coupling constants, etc. In the standard approach they are dependent on the cut-off 
momentum and it is tacitly understood that the experimental dependence of these parameters 
on the squared transferred momentum Q 2 is equivalent to the dependence on the cut-off. In 
our wavelet-based approach the g = g(a) may be the coupling constant of the given field <j) a (b) 
of the given scale a, rather than the effective coupling constant for the harmonics of all scales 
up to a. For this reason different notations are used hereafter in the Fourier and the wavelet 
representations of field theories. 

2 Wavelets and scale-dependent functions 

2.1 What is wavelet transform 

Wavelet transform entered mathematical physics from geophysical applications J31 H] as a prefe- 
rable alternative to the Fourier transform in the case when localization in both position and mo- 
mentum is simultaneously required. Formally the Fourier transform and the wavelet transform 
are on the same footing: both are decomposition of a function with respect to the representa- 
tions of a Lie group. However, the former uses Abelian group of translations, the latter uses 
non-Abelian group of affine transformations 

x' = aR(6)x + b, R(9)£SO d , x,x',b£R d . 

Let us remind the general formalism. For a locally compact Lie group G acting transitively on the 
Hilbert space 7i it is possible to decompose vectors <f> G 7i with respect to the square-integrablc 
representations U(g) of the group G 0IEJ: 



C^ / \U(jg)il>)dM(1>U*(jg)\<l>), 

JgeG 
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where d[i(g) is the left-invariant Haar measure on G. The normalization constant C^ is deter- 
mined by the norm of the action of U(g) on the fiducial vector tp € TC, i.e. any ip £ 7i that 
satisfies the admissibility condition 

C i , = U\\- 2 f mU{g)^)\ 2 d^{g)<^, 

JgeG 

can be used as a basis of wavelet decomposition (J2J. 

Hereafter we assume the basic wavelet tp is invariant under SOd rotations ip(x) = ip(\x\) and 
drop the angular part of the measure for simplicity (R(0) = 1). After this simplifying assump- 
tion, the left-invariant Haar measure on affine group is dfj,(a,b) = a d+1 . The representation 
U(g) induced by a basic wavelet ip{x) is 

g:x' = ax + b, U{g)ij){x) = a' d/2 ip ( ^— ^). (3) 



Therefore, in the Hilbert space of square- integrable functions L 2 (M ), with the scalar product 
(/>#) := / f( x )g( x )dx, a function <fi G L 2 (R d ) can be decomposed with respect to the represen- 
tations of affine group © : 

Mx) = c -,j a -^{^ Ub) ^, (4) 

The coefficients of this decomposition are 

Mb) = f a- d ' 2 i> (—) 4>{x)d d x. (5) 



The d-dimensional bold vector notations are dropped hereafter where it does not lead to a con- 
fusion, and the basic wavelet ip{x) is assumed to be isotropic. Here we use normalization for the 
rotationally invariant wavelets 



S d \k\ d 
where the area of the unit sphere in d dimensions Sd has come from rotation symmetry. 

2.2 Scale-dependent functions 

Up to this point we considered only a projection of a (j) G L 2 (M ) functions onto the basis con- 
structed of shifts and dilations of the basic wavelet ip(x). If we substitute the inverse wavelet 
transform © into any physical theory we have started with we will reproduce it identically. 
However, the wavelet coefficients <p a {b) may have their own physical meaning. The convolu- 
tion © can be considered as a "microscope" that scrutinizes a field or a signal 4>(x) at a point b 
at different resolutions. In this sense the function ip can be considered as an apparatus function 
of the measuring device. 

The most known case of such theory is the measurement of the turbulent velocity field: the 
mean velocity and the PDF of fluctuations of different scales should not be the same ^ [7j. 
Experimentally averaging of the field in a physical volume L is usually described as 



v L = L~ a / K(x - y)v(y)d d y, 

JL d 

where K is some averaging kernel. But the "no-scale" field v(y) may not exist at all: consider 
a mean velocity of fluctuations at a scale L less than a mean free path of the particle. 
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That is why the wavelet coefficients (p a (b) may have their own operational meaning, even if 
the "no-scale" function (j)(x) does not exist. In this case a sum of all fluctuations with scales 
equal or greater than a given scale L can be defined [S]: 



2 f°° wo fx-b\ ,.dad d b 



**>&-§; F'-"*^)*-® 



a d+i 



This is a close analogy to the Wilson's RG approach to critical phenomena |9j I10| ITT], where 
a magnetization (or velocity, or other order parameter) of scale L is taken in the form l|Tjl. 

2.3 Wavelets and random processes 

If a function to be analyzed by continuous wavelet is a random function, its wavelet transform 
is also a random function. So, instead of the usual space of the random functions f(x, •) € 
(Q,A, P), where f{x,ui) £ L 2 (IR rf ) for each given realization uj of the random process, we can go 
to the multiscale representation provided by the continuous wavelet transform (|5]l: 



The inverse wavelet transform 

reconstructs the common random process as a sum of its scale components, i.e. projections onto 
different resolution spaces. 

The use of the scale components instead of the original stochastic process provides extra 
analytical flexibility of the method: there exist more than one set of random functions W(a, b, •) 
the images of which have coinciding correlation functions in the space of f(x, •). It is easy to 
check that the random process generated by wavelet coefficients having in (a, k) space the 
correlation function 

(W(ai, h)W(a 2 , k 2 )) = C v ,(27r) d ^(A;i + k 2 )a d+1 5(ai - a 2 )D 

has the same correlation function as white noise has: 

(f( Xl )f(x 2 )} = D^ d { Xl - x 2 ), 
(f(h)f(k 2 )) = (27T) d D S d (k 1 + k 2 ), 



(W V ,(ai,A;i)W V) (a 2 ,A; 2 )) = {2^) d D 5 d {k l + M(aia 2 ) d / 2 ^(ai&i)<Ka 2 & 2 ). 

Therefore, the space of scale-dependent random functions is just richer than the ordinary space 
of random functions: since we have an extra scale argument a here, we can play with the PDF 
and correlations of random functions <j) a (b, ■) to achieve required properties in ordinary space, 
holding at the same time some other limitations, say, on singular behavior. 

3 Stochastic quantization 

3.1 Parisi-Wu stochastic quantization 

Quantum field theory, as is understood in technical sense, is a tool to calculate vacuum ex- 
pectation values (v.e.v.) of the field operator products {4>{x\) . . . 4>(x n )) - the Green functions. 
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Usually it is done in the functional integral formalism: if the action of the field 4> S[(j>] is known, 
the v.e.v. can be derived by taking functional derivatives of the generating functional 

W[J(x)]= [ V4>e lSm+l $ J{x) ^ x)dx (6) 



with respect to the formal source J{x). Unfortunately, infinite perturbation series obtained by 
such differentiation often give infinite results and require regularization. 

However, there are alternatives to this method. Since any quantum system interacts with 
environment, it can not be in a pure state, and averaging over all field configurations should be 
performed taking into account averaging over the states of environment. This can be done by 
considering the system environment as a thermostat, described by a Gaussian random noise; the 
desired v.e.v. state will be a result of relaxation to the thermodynamic limit. Thus, instead of 
the Feynman path integral ©, one can operate with stochastic differential equations. Because 
of equivalence between Euclidean QFT and the stochastic problems, the stochastic methods 
have found wide applications in both numeric and analytic solution of QFT problems, especially 
in regularization. Among those, the method of stochastic quantization of gauge fields, proposed 
by G. Parisi and Y. Wu |12| . have been attracting attention for more than 20 years. The idea of 
the method is as follows. Let <Se[^>] be the action Euclidean field theory in Mr. Instead of direct 
calculation of the Green functions from the generation functional of the field theory, it is possible 
to introduce a fictitious time r, make the quantum fields into stochastic fields (f>(x) — > 4>(x,t), 
x G M d , r£l and evaluate the moments (0(xi,ri) . . . 4>(x m ,T m )) by averaging over a random 
process <f>(x,T,-) governed by the Langevin equation 

X - T) " S «X,T). (7) 



dr 84>(x,t) 

The Gaussian random force is 5-correlated in both the M d coordinate and the fictitious time r: 

(rj(x, r)rj(x', r')> = 2D S(x - x')5{t - t'), ( v (x, t)> = 0. (8) 

The physical Green functions are obtained by taking the steady state limit: 

G(xi,...,x m ) = lim ((j)(xi,T)...(j)(xm,r)). 

T— >0O 

The stochastic quantization method is much preferable to ordinary methods for it does not 
have problems with gauge fixing and does not require incorporation of higher derivatives as 
continuous regularization methods do. 

Concerning the renormalization of stochastically quantized theory, the matter practically 
comes out to the renormalization of the Langevin equation Q in (d + 1) dimensions instead 
of renormalization of the original theory with the action functional S[(p] in d dimensions. The 
renormalization of the Langevin equation is usually done by constructing the characteristic 
functional 



Z H = («p(/*W(„ TW („r: 



solutions 



where the statistical averaging (• • •) is taken over all solutions of the Langevin equation l|7jl. 
The summation over all solutions in functional integral formalism is achieved by expressing the 
functional 5-function as a functional integral over an imaginary auxiliary field (p: 

Jd<t>(x,T) 5S . A /" j j/00(x,t) SS , A 
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This is so-called Martin-Siggia-Rose field-doubling formalism |13] . The resulting theory has the 
form of a field theory with constraints in (d + 1) dimensions. Under the change of the cut-off 
scale all fields, coordinates and parameters in such a theory are transformed according to their 
canonical dimensions. 

The renormalization of such a theory has much in common with the renormalization of gauge 
theories: for the constraint © is similar to gauge fixing, and the functional determinant arising 
from the this term under the change of variables results in ghost fields and BRST symmetry. 
The details of renormalization of stochastically quantized theories can be found e.g. in |14| ll5j. 

3.2 Wavelet-based stochastic quantization 

Taking into account the above mentioned flexibility of stochastic processes defined in the space 
of wavelet coefficients, it appears to be attractive to use wavelet-defined noise for stochastic 
quantization. Applying the wavelet transform (in spatial coordinates, but not in fictitious time) 



-d' ' x 



4> a {b,-) = J a~^[-—U{x,-)d a x 

to the fields and the random force in the Langevin equation Q, we get the possibility to change 
the white noise (JSJ) into a scale-dependent random force 

(fjaAkun^ih,^)) = C^) d 5 d {k x + k 2 )5(n - T 2 )a x 5{ ai - a 2 )D(a 1 ,k 1 ). (10) 

(Here and after we use CWT in L 1 norm instead of L 2 .) 

In the case when the spectral density of the random force is a constant D(a\, k\) = Dq, the 
inverse wavelet transform 

<t>( x ) = TT — I 77r^7^exp(z(/cx-u;T))^(afc)c/> a (A;,u;) (11) 

Cy> Jo a J ( 27r ) 27r 

drives the process (jlUj) into the white noise (|5|). 

In the case of arbitrary functions <j> a (x, •) we have more possibilities. In particular, we can 
define a narrow band forcing that acts at a single scale 

D(a, k) = a 5(a - a )D . (12) 

The contribution of the scales with the wave vectors apart from the the typical scale a^ is 
suppressed by rapidly vanishing wings of the compactly supported wavelet ij){k). 

Here we present two examples of the divergence free stochastic perturbation expansion: (j) the 
scalar field theory (p s , (ii) the non-Abelian gauge field theory. 

3.3 Scalar field theory 

Let us turn to the stochastic quantization of the ^-theory with the scale-dependent noise |16j . 
The Euclidean action of the c/) 3 -theory is: 



S E [<t>(x)} = / d d x 



iW)' + V + ^ 



The corresponding Langevin equation is written as 

9HX ' T) + ( ~ A^ + m^ + *A = n(x, r). (13) 



Or V 2! 
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Substituting the scale components in representation (|llj) into the Langevin equation ((T3]) we get 
the integral equation for the stochastic fields: 

. , 2 2n / /i \ /i \ A ~ . — / 2 \ f d d k\ dijj\da\dai 

(—lui + k +m )<pa{k, to) = 7] a {k, to) - —ifr{ak) — 



2 y VCv/ J (2vr) d 2vr a\ a 2 ' 
ij){aiki)i>{a2{k - ki))(fi ai (ki, ojx) (pa 2 (k - h,u -u>i). (14) 

Starting from the zero-th order approximation (fro = Gorj, with the bare Green function Go(k, uj) = 
l/(—iuj + k 2 -\- m ), and iterating the integral equation (|T4*|) . we get the one- loop correction to 
the stochastic Green function: 

G{k,u) = G (k,uj) 

+ \ 2 G 2 (k,u) f ^L^-2A(q)\G (q, n)\ 2 G (k - q,u - fi) + • • • , (15) 

where A (A;) is the scale- averaged effective force correlator 

A(fc) = ^-/ -\ifr(ak)\ 2 D(a,k). 
L"ip Jo a 

In the same way other stochastic moments are evaluated. Thus, the common stochastic diagram 
technique is reproduced, but with the scale-dependent random force (|10j) instead of the standard 
one (jSJ). The diagrams corresponding to the stochastic Green function decomposition ((To]) are 
shown in Fig. ^ 

Q 

K K K J#^ K 

—h- ■ — > — K3^ + 

Figure 1. Diagram expansion of the stochastic Green function in </) 3 -modcl. 

It can be seen that for a single-band forcing Q12JI and a suitably chosen wavelet the loop 
divergences are suppressed. For instance, use of the Mexican hat wavelet 

$(k) = {2it) d l 2 {-ik) 2 exp(-k 2 /2), Gp = (27r) d , 

and the single-band random force (|12jl gives the effective force correlator 

A(g) = (a g) 4 e-(°°9) 2 A,. (16) 

The loop integrals, taken with this effective force (jlfij) . can be easily seen to be free of ultraviolet 
divergences at q — > oo: 

G 2 (k,u;) = G (k,u;)J j^y d 2A(q) J^ — - 2 - ^ - -^ -^ _ - - (jfe _ -^ - - 2 . 

3.4 Non-Abelian gauge theory 

The Euclidean action of a non-Abelian gauge field is given by 



S[A] = ±Jd d xF« l/ (x)F^(x), 



F; u {x) = d,A a v {x) - d v A a Jx) + gt bc A b Jx)Al( X ), (17) 
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where / are the structure constants of the gauge group, g is the coupling constant. The 
Langevin equation for the stochastic quantization of gauge theory (|17|) can be written as 

dA a (x t) 

/1 g r ' + {-^ud 2 + d^d v )A a u {x,r) = <(x,r) + t£(x,r), (18) 

where rftix, r) is the random force and U^(x,t) is the nonlinear interaction term 

The two terms standing in the free field Green function correspond to the transversal and the 
longtitudal mode propagation: 

Similarly to the scalar field theory, we can use the scale-dependent forcing (|19fl in the 
Langevin equation (|18[). Since there is no dynamic evolution for the longtitudal modes in the 
Langevin equation 1)18(1 . it is natural to use the transversal scale-dependent random force 

«#l,n)77L 2 (fc2,r 2 )) = (27r) d ^(fci+fc 2 )5(r 1 -r 2 )T^(fc 1 )^a 1( 5(o 1 - a 2 )D(a 1 ,k l ). (19) 

Let us consider a gluon loop with two cubic vertices. Summing up over the gauge group indices 
(§5) f abc dbdf der $cr = if<5 ae C2, with C 2 = N for SLOv gauge groups, we can write the gluon 
loop as a sum of two diagrams - those with the transversal and the longtitudal stochastic Green 
functions: 

//70 d d n 
— j^- I N I (k,u J , q ,n)lj lu (k, q )2A( q ), 
1=1,1, n *> ' 

where 

1 



N(k,q) 



-zQ + q 2 






l(iu(k,q) = V^K\(k, k - q, q)T Xl {q)V avl {k - q, k, -q) [ L ( k _l 



As it can be observed after explicit evaluation of the tensor structures l^ v and lj^ u using the force 
correlator 1)19(1 . and integration over d£l, the wavelet factor in the effective force correlator A(q) 
suppresses the divergences in the case of a narrow-band forcing (|12j) . The power factor k n of the 
basic wavelet ip that provides ^(0) = and makes the IR behavior softer. In this respect the 
wavelet regularization is different from the continuous regularization f d yRh(d 2 )r](y,T), which 

_fe 2 __ 

makes UV behavior softer by the factor e a 7 , but do not affect the IR behavior, see e.g. [17j . 

4 Wavelet-based Euclidean field theory 

4.1 QFT on a Lie group 

Let us consider a Euclidean field theory determined by a characteristic functional 



W E [J] =M fvcj)exp 



Se[<P(x)] + / d d xJ(x)cj>(x) 
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where Se[ a ( x )] is Euclidean action, N is formal normalization constant. The (connected) Green 
functions (m-point cumulative moments) are evaluated as functional derivatives of the logarithm 
of generating functional W.e[J] = e~ E ^ ': 



G m (x 1} . . . , X m ) = (<f>(xi) . . . 4>{Xm)) 



5J(xi) ■■■8J(x m ) 



\uW E [J}. 



j=o 



Applying a formal partition of a unity |2j) with respect to a Lie group G we yield a theory with 
the generating functional 



W g [J{g)]=N Jv<f>(g)exp 



S g [<f>(g)] + / dn{g)J{g)4>{g) 



G 



g^G 



and appropriately defined Green functions. 



4.2 Wavelet-based action 

Let us consider the theory of a massive scalar field with polynomial interaction 



S E [ 



d d x 



l^f + ^ + ww 



(20) 



that can be alternatively interpreted as a theory of classical fluctuating field with the Wiener 
probability measure VP = e~ SE ^T>4>. In this case mr = \T — T c \ is the deviation from critical 
temperature and A is the fluctuation interaction strength. 

To introduce the scale-dependent fields <p a (b) into the action we can just substitute the 
fields <j) into (|20j) using wavelet transform (|lj). For convenience and recording purposes, we 
rewrite continuous wavelet transform (JSJ), © using the L 1 norm 



Mb) 

<f>(x): 



1 , fx-b 



or \ a 

2 r°° da 

C,i, Vn a 



(x)d d x, 



'i> Jo 



d d bMb)- d i>(— 
ar \ a 



(21) 



C i> - S d 



&*k-> 



|V>(«fe)l : 



da < oo. 



This provides the fields <j)(x) = (x\(f>) and <p a {x) = {a,x;ip\(f)) with the same physical dimension. 
The generating functional for the scale-dependent functions 4> a (b) is 



w[u 



VUb) ( SEwlMb)} + I Ja(b)Mb)^-) 



where J a (b) is a formal source ("external force"), which corresponds to the fluctuations of given 
scale a localized near a given point b. The corresponding Green functions a given by 



G(ai,xi;...;a n ,x n ) 



S\aW[J a (x)] 



SJa 1 {xi)5Ja n {x r 



y «i 



(x 1 )...(p an (x n )). 



J=0 



Let us consider a theory with polynomial interaction V((f>) = —^ 



3. An 



Sew 



Tfl 

-(t> ai (bi)D(a 1 ,a 2 ;h -b 2 ) A a 2 (h) + ~^-\ A a( 



9 



+ -.V, 



n 



ra\...a r 
I "6i...6n 



(t> ai {bi)---(t>a n (b n ), 



(22) 
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KX" = J d d *V( Xl , . . . , x n ) JI ±* (^) , (23) 

where integration over all pairs of matching indices -pi — *■ is assumed. The inverse propagator 
is the wavelet image of (— d 2 + to 2 ), that is 



D(a±, a,2, k) = i/j(aik)(k + to )ip(a,2k). 

As it is seen from (|23j) . the local interaction term (j) n becomes nonlocal after the application of 
wavelet transform (|21|). Of course, since the wavelet transform provides the partition of a unity, 
the integration over all scale arguments in the theory (|2'2j) drives it back to the original field 
theory with polynomial interactions, which is UV divergent. However, using the Wilson's RG 
ideas we shall show how the nonlocal wavelet theory (1221) can be made into a local one, which 
has UV finite behavior. 

Let us suppose that we manage to make the wavelet theory with nonlocal interaction l)23|) 
into a local one with the coupling constant explicitly dependent on scale g = g{a). The simplest 
case of the fourth power interaction of this type is 

T . r ,i 1 I q(o>) ,A, r ,dadb , . . ,, 

Vint[<t>] = -^j -Y<t>a(b) — , g{a) ~ a". (24) 

The one-loop order contribution to the Green function G2 in the theory with local interac- 
tion (|24|) can be easily evaluated j!8j by integration over a scalar variable z = ok: 



I 



a u \i>(ak)\ 2 d d k da _ (l/) f d d k k~ 



k 2 + to 2 (2ir) d a * J (2vr) d k 2 + to 2 

where 



- 1 Jm 



C^ = C^, 1 I Mz)\ 2 z l 



Therefore, there are no UV divergences for v > d—2. However, the positive values of v mean that 
the interaction strengths at large scales and diminishes at small. (This is a kind of asymptotically 
free theory that is hardly appropriate say to magnetic systems.) 

We have to admit that the idea of using wavelets for regularization and blocking of degrees 
of freedom has been considered by G. Battle ^1 an d P. Federbush [20] in lattice settings. The 
idea was to sum up the discrete wavelet expansion of a function starting from a large infra-red 
scale (set to unity) up to the infinitesimally small scale instead of integrating in Fourier space. 
The expansion was taken in the form 

4>{x) = ^akUkix), u k (x) = ^ k(x), (25) 



v^AT 



TO 



where k = (n, m, s) is a multiindex that incorporates translations and binary dilations on the 
lattice as well as internal degrees of freedom of the field <j), ^>k( x ) forms an orthogonal basis 
in L 2 (M ). It have been understood then that different subgrids involved in the wavelet expan- 
sion (|25|) represent different independent degrees of freedom [^J, rather than being just different 
approximation of the same field like that in Kadanoff block expansion. The corresponding dia- 
gram technique, called "wavelet diagrams" and originally proposed by Federbush, can be found 
in Battle's book [15] . 

For a particular type of wavelet basis, the Lemarie wavelets, it was shown that due to the 
big number of vanishing momenta of the basic wavelets ^ the correlations between different 
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block variables {akOt-k') decrease as a high inverse power of the separation between blocks. As 
a consequence of this a field theory (S , [0(x)] becomes UV-finite without any renormalization. 
However the decoupling of interaction between different scales seems to have quite general 
nature and that is the reason we exploit this idea in present paper in terms of continuous 
wavelet transform without using any lattices. 

Now we will try to show how the local interaction of type l)24j) can emerge. Following 
K. Wilson ^3 ^J we consider a system of spins with spacing Lq. The idea of RG applied to 
a spin system is that the interaction of spins separated by a distance Lq, can be substituted by 
interaction of blocks (of 2 d spins in each) of size 2Lq with the high-frequency details absorbed 
into the interaction constants. The same procedure is then applied to the blocks of 2Lq size, 
blocks of 4Lo size etc. - this is so-called Kadanoff blocking procedure. On each step of the 
blocking procedure the transformation from the scale L to 2L is done by integrating over the high 
frequency degrees of freedom K > -j- with appropriate adjustment of the coupling constants. 

In continuous limit, i.e. in Landau theory of ferromagnetism, when spins are averaged into 
the mean magnetization <p(x), the energy of the spin system is given by the Ginzburg-Landau 
free energy functional 



F E [4>]= I d d x 



^) 2 + f0 2 + ^> 4 



(26) 



where integration is performed over the domain occupied by the spin system. The field theory 
(|26l) gives divergences in correlation functions. To cope with it we remind that <j>(x) is the 
mean magnetization of a macroscopically big block, i.e. one that contains enough spins to make 
statistical averaging valid. This means that the field 4> should be substituted by one containing 
only the fluctuations of size larger than a typical scale L: 

So, the original Ginzburg-Landau functional is made into effective action for the large-scale 
fields 0rw 



Fl[4>( L) ] = fd d x V-{d^ L) ? + ^<pf L) + ^4>f L) 



(28) 



in which the effects of small-scale fluctuations with \k\ > 2£ are absorbed into the coupling 
constants Rl and Ul- It should be noted that the projection (|27|). performed by filtering out 
the Fourier harmonics with \k\ > -j-, is not the only possible type of smoothing: an exponential 
or proper-time cutoff can be used as well [22| . Similar projection can be obtained if we apply 
the wavelet transform with certain given kernel <p{x) — > 4> a (b)- 

Wilson suggested an elegant way to determine the dependence of Rl and Ul on the cutoff 
scale L by averaging over fluctuations in the shell [L,L + 6L) ^OlEH^D- Since we understand 
the fluctuations of a given scale 4> a {°) ~ now in wavelet notation - as physical fields measured 
at scale a, we will reproduce the Wilson's derivation integrating the scale-dependent free-energy 
F[(f) a (b)] in appropriate limits over the logarithm of scale a~ 1 da. The larger is the scale L, the 
more fluctuations of smaller scales contribute to free energy. 

Let Lq be the smallest size of the system - the distance between spins in the theory of 
ferromagnetism. Then, the Kadanoff blocking procedure is a mapping 

H^(R , Uq) - H^(R U Ui) - HM(R 2 , U 2 ) - • • • , 
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where {ip k } is the basis for the finest resolution Hilbert space, with H^°'(Rq,Uq) being the 

Hamiltonian acting on this space, {tp k } is the basis of the next coarse-grained space of the 
resolution 2Lq, etc. At each step of the coarse graining process some details are lost, so any 
function ip k ' of the basis {tp k } can be expanded in the basis {ip k }; and similarly for all next 
scales: 



^ = Ei + M 



since {ip k } is more detailed than {ip k }. 

Expanding the more coarse-grained basic functions in terms of the less coarse-grained one we 
obtain the expansion of the parameters of the coarser Hamiltonian H^ +1 > in terms of the finer 
Hamiltonian H^'. Wilson suggested a qualitative way of doing that. Suppose we know the free 
energy functional Fl [4>] of scale L and need to calculate the next coarser functional Fl + $l [<ft] of 
scale L + 5L. Then, we have to add fluctuations of all scales within the range [L, L + 5L) to the 
original theory and integrate over those fluctuations. If ift(x) is a normalized basic function of 
scale L, e.g. a wave packet with k ~ 2j£, such that 



<fx\^(x)\ z = 1, / d a x\dijj(x)\ z ~ L~ z , (29) 

we just make 4>{x) — > (j)(x) + ap(x) and integrate over the scalar amplitude c: 

/oo 
dme- F ^ +m ^. (30) 

-oo 

Substituting the free energy (|28|) into the condition (|3U|) and taking into account the condi- 
tions (|29|) we get 

I™ dme-^ +R ^ U ^) . (31) 

J — oo 



e -F L+SL [4>] =e -F L [ 



The derivation is not very rigorous, and referring the reader to the original paper jllj for the 
details, we just have to state that using the formal properties of Gaussian integration in (|3*Tj) . 
up to the numerical factors, we get a logarithmic contribution to the free energy 

Fl+sM = F L [<f>] + l - In f-1 + R L + \ul4> 2 ^ ■ (32) 

To derive the dependence of Rl and Ul on scale L two more tricks were used by Wilson: 
(i) the logarithm in 1)32(1 is decomposed into a Taylor series; (ii) integration over the phase space 
volume corresponding to the d.o.f. related to the wavepacket ip is changed into integration over 
the coordinate volume 5V using the uncertainty principle. 

The logarithm is decomposed into the power series only up to the second order - to get the 
forth order in the fields. Up to the terms that do not depend on cj) we get 

l - In Q^ + R L + \u L (jA = -^ (u L L 2 <p 2 - ^LU 4 ~ RlUlLU 2 ) +■■■■ (33) 

The next trick is to count the degrees of freedom corresponding to the localized wavepacket ip. 
Since the scales are in the range [L, L + SL), the volume in momentum space is ^ ~ 77 • The 
phase space volume, since p ~ 27r/L, is 

5T ~ Vp d - 1 Sp ~ VL~ l ~ d 5L, 



Scale-Dependent Functions, Stochastic Quantization and Renormalization 13 

and, hence, the coordinate volume per one degree of freedom is 

L i+d 



SV 



SL 



Thus the equation IJ33JI can be formally multiplied by 1 = SVL 1 SL and integrated over the 
volume SV. This results in equations 

R L+ 8L = Rl + \ (L l - d U L - L^ d R L U L ) SL, 
U L+ s L = U L - h 3 ~ d UlSL. 

See Uni for general theory of renormalization in field theory and critical phenomena. 

The above Wilson's consideration was presented only to show how integration in the thin 
shell below the cutoff renormalizes the parameters of the action. More rigorous formalism that 
accounts for the change of the action parameters by means of integration over the infinitesimally 
thin shell 1 — SI < \q\ < 1 in momentum space, followed by rescaling, consists in constructing 
differential equations for the variation of the action functional. It is known as exact renorma- 
lization group (ERG). The ERG equations are rather complicated and can be found elsewhere 
|24| I10| 125] . If solved numerically, the ERG equations could provide exact scale dependence 
of the coupling constants on cutoff scale in terms of the considered model. The problems are, 
however, that (i) it is often impossible to solve these equations, (ii) the cutoff scale is not the 
same as the scale of observation, and hence it is not obvious how to interpret the ERG results; 
(hi) and the last, but not the least, that it is not clear how to separate the modes to be integrated 
over from those to be renormalized by this integration using the space of functions that depend 
only on momentum (p(k) (JO], Fig. 11.2). Here we propose an alternative solution of the mode 
separation problem: we just extend the space of functions adding the scale argument explicitly 
#■)-&(■)■ 

Having the information on the dependence on scale of the effective mass Rl and coupling 
constant Ul we can address the question, how the effective theory Q28J1 that comprises nonlocal 
interactions of all modes larger than L can be transformed into a theory with the local interac- 
tions of scales of the type (|24|). To do this we approximate the original Ginzburg-Landau action 
F, by a new action 



F E [ 



dad d b 



\4> ai {b{)D<i> a2 {b 2 ) + ^\Ub)\ 2 + g ^§- 



(34) 



That meets physical requirements by construction: for a free field (g(a) = 0) it coincides with 
the original free action; for the interacting fields, it just satisfies the Wilson's assumption that 
only the fluctuations of the close scale interact: in the case of the action (|34[) . certainly the 
equal scales only. 

Since the action (]34j) explicitly involves integration in scale variable da/a, we can attribute 
the difference in the free energies 1)33]) to interaction of the appropriate "renormalized" modes, 
governed by the equation (|34|). in the shell [L, L + SL). 

The n > 2 the polynomial interactions are nonlocal in wavelet space. To derive their interac- 
tion constant dependence on scale, we use the physical assumptions that only the fluctuations 
of close scales can directly interact to each other. So, in the self-interaction of the field <f> 



U L I <f>(x) n d d X = U L I Vg^aM ■ ■ ■ ^ an (bn) f[ 
J Lj J Li ■ -i 



dajd d b; 
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the terms with | lnaj — lnaj| <C 1 should dominate, from where we can postulate that 

Ul f V£;;£K(h)---ct>aMfl^^= f - f g(a)\Mb)\ n d d b (35) 

JL ■ -, a i J L a J 

for some g(a). 

It should be noted that similar idea of orthogonalization of fluctuations of different scales has 
been already applied to Ginzburg-Landau model by C.Best using discrete wavelet transform. 
However, the Best's paper |26j . having started from continuous Ginzburg-Landau model, uses 
orthogonal wavelets and sufficiently strong assumptions that fluctuations of different scales are 
^-correlated in both position and scale (equation (4) from J2H])- This is definitely not the case for 
arbitrary (non-orthogonal) wavelets and standard assumption of Gaussian nature of fluctuations 
in ordinary continuous model (|26|). Nevertheless, the Best's idea is numerically convenient and 
it is interesting for what type of wavelets such behavior is really observed. 

Let us perform the asymptotic estimation of the behavior of g(a) at the known behavior 
of Ul- Let us assume that interaction takes place at the thin shell of scales [L, L + SL) and set 
a\ = ■ ■ ■ = a n = L. Since we need the behavior of the local interaction of the type (|24|). we set 
&i = • • • = b n = in the r.h.s. of (|33|) to get the upper bound for g(a). Doing so, we get 



tft.o* = Jd d x^r (~) = a^ d Jd d yr(y). 



In Wilson's theory (n = 4, d < 4) 

Ul = YJi^d ~ L 

2^T + const 

and hence 

g{L) - L~ 3d U L - L~ 2d - 4 . (36) 

The important point is that in the theory (|34|) the observable quantities are the correlators of the 
scale-dependent functions {4> ai (bi)(j) a2 (b2))- If we deal with "renormalized" action (|34|) . where 
only equal scales interact polynomially, the loop integrals can be easily evaluated for a fixed 
value of scale. In the fourth power interaction theory the one-loop contribution to the Green 
function will be 

da , . f d d k \4>(ak)\ 2 f a ™* da ,.,.,. 

t 9(o) J w?ikh = jL t 9(o)l(o) ' (37) 

The ^-dependent integral in the last equation (|37|) can be explicitly evaluated. As an example, 
let us present a "bad" case of the Morlet wavelet ip(k) = exp ( ikz — |- J , which is not admissible 
in the sense C^ = oo, 

rco e -aV gd-i d d d=4m=1 - a * e °?Ei(l, a 2 ) + a 2 a 2 ~ 



2 



q 2 + m 2 (27r) rf a 

For admissible wavelets (C^, < oo) the IR behavior is milder. 

Certainly, the integral J Qmax —g{a)L(a) in any positive range of scales and with the coupling 
constant g(a) taken from (|36(i is finite. 
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5 Conclusion 

In this paper we attempt to construct a theory of scale-dependent fields 4> a {b) starting from an 
action functional that explicitly depends of these fields, rather than being a local functional S[(j)] 
truncated at certain scale L by averaging over the short-wave fluctuations \k\ > -j£. The diffe- 
rence is that in our theory it is meaningful to consider the correlation functions {cj) ai (oi)0a 2 (&2)) 
between different points and different scales. This is important because any experimental mea- 
surement is performed not exactly in a point, but in a certain vicinity of a point, the size of 
which is constrained at least by the uncertainty principle: the higher momentum is used in the 
measurement, the smaller is the vicinity and the stronger is the perturbation caused by the 
measurement. Thus we have a strong need do construct a theory in terms of what is really mea- 
sured in experiment - scale-dependent (wavelet) fields 4> a (b), rather than in terms of abstract 
"no-scale" functions (f>(x). 

Standard RG approach meets this need only partially: it enables to study the dependence of 
correlations on separation between points, b = b\ — b% in our terms, but not on the typical scale 
of interaction. The latter is described in terms of the correlation length £ = £(T), which has 
a universal behavior near the critical temperature £(T) ~ (T — T c )~ u . The Kadanoff hypothesis 
that blocks of spins interact exactly as spins themselves becomes strictly valid only in critical 
regime £ — ► oo at T — ► T c . In other cases there are no reasons to assume that fluctuations of 
different sizes do not interact to each other and (4>a 1 (bi)4 , a 2 (b2)) = for a\ ^ ai, instead the 
inter-scale interaction is expressed by wavelet diagrams. 

As it concerns the criticality, our approach may be close to the ideas of functional self- 
similarity (FS) [23 : at the critical point the blocks of spins can interact with each other in 
the same way as primitive spins themselves, but the coupling constant is scale-covariant, i.e. its 
transformation with changing of scale depends on the value of coupling constant, but not on the 
scale. 

For a quantum field theory models, which arise from the high energy physics problems, the 
Kadanoff block-averaging procedure can not be directly applied. Instead the universal behavior 
of the coupling constant g = g(a) with respect to the scale transformations a' = e x a should be 
considered as a fundamental symmetry of the physical system. Indeed, the coupling constant 
g(a) is a charge of a given particle with respect to a given interaction measured at a given scale a. 
Since the experiments can be performed at different scales ao < «i < «2 < • ■ ■ , we should expect 
a universal relation between the values of charges at different scales, and it should be an invariant 
charge, such that g = C(ao,eo) = C(o2/ao,e 2 ) = C(oi/oo,ei), or, since e\ = £(02/01,62), 

C(a,g) = aa/t,C(t,g)). (38) 

The transformation of the charges between different scales (|H8jl has a structure of the Abelian 
group Ttjt 2 = 2fj o Tt 2 . It describes the evolution of the coupling constant with the change of 
measurement parameter - the resolution a. 

In differential form the self-similarity equation IjMHjl is most naturally expressed in the loga- 
rithmic coordinate I = log a: 

G(l + X,g) = G(l,G(X,g)), 



or 



9 /%)#Wa)=0, where (3(g) - dG{X > 9) 



dl dg J d\ 

The latter is referred to as the Lie-Ovsyannikov equation. 



A=o 
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In general case a quantity f(l,g), which is not an invariant charge, is not nullified by the 
generator of RG transform e~^ R f(l,g) 7^ f{l,g), but the generator 

determines the evolution of physical quantities under the scale transformations. 

In our approach we extend the space of functions <j){x) E L 2 (M d ) to the space of multiscale 
(wavelet) fields 4> a (b) for which both the position b and the scale a are considered as independent 
variables. The RG transform therefore becomes just an Abelian subgroup of the whole symmetry 
group of the action S[<p a (b)] described above in this paper. The difference between our approach 
and the exact renormalization group approach (ERG) consists in the fact that ERG generating 
functional 



Z A [J] = / II ^(P)e- S[<M1 - /<W 



|p|<A 

is just a sum of all fluctuations truncated by a cutoff momentum A in their Fourier repre- 
sentation, and therefore having certain symmetries related to the rescaling of A - the cutoff 
parameter. In our approach the logarithm of scale I = log a becomes an independent variable, so 
the technique of Lie symmetry analysis can be applied with respect to I on the same footing as 
with respect to the position x. In this way the vector field that generates arbitrary transforma- 
tions of the action S[4> a (b)] due to the transformation of dependent and independent variables 
takes the form 

( d d\ d d 

We hope that such symmetries could be used to find solutions of the wavelet renormalization 
equation (|35|). 
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